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Abstract 

We develop linear discretization of complex analysis, originally introduced by R. Isaacs, 
J. Ferrand, R. Duffin, and C. Mercat. We prove convergence of discrete period matrices 
and discrete Abelian integrals to their continuous counterparts. We also prove a discrete 
counterpart of the Riemann-Roch theorem. The proofs use energy estimates inspired by 
electrical networks. 
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1 Introduction 

The history of linear discretizations of harmonic and holomorphic functions can be traced back to 
the early finite-element literature [8, 10]. The discretization of the Dirichlet variational principle 
leads to linear discretizations of the Laplace equation and of the Cauchy-Riemann equations on a 
^ . square grid. Lelong- Ferrand developed the theory of discrete harmonic and holomorphic functions 
7— ( \ on a square grid to a level that allows her to prove the Riemann mapping theorem [13, 16]. Duffin 
^ [ takes first steps to extend the theory beyond the square grid. He considers arbitrary triangulations 
O ' °f planar domains and discovers the now famous cotan- weights [11]. The cotan-Laplace-operator 
q ■ turned out to be important in discrete differential geometry of surfaces [ ] and in computer 
graphics [21]. 

A striking feature of some 2-dimensional discrete models in statistical physics is that they 
exhibit conformally invariant properties in the continuous limit: for site percolation on a triangular 
grid [28], for the random cluster model [: ], for the Ising model [6], for domino tiling [15]. In all 
cases, the linear theory of discrete analytic functions on regular grids was instrumental. 

Mercat generalized the linear theory from planar domains to discrete Riemann surfaces [17] and 
in particular introduced discrete period matrices [19]. Numerical experiments to compute period 
matrices for triangulated surfaces and compare them with known period matrices for the smooth 
surfaces were performed in [3, 2]. First attempts to prove the convergence were made in [18] but 
it remained an open problem; see Remark 6.5. 

A different "nonlocal" linear theory for discrete holomorphic functions on triangulated surfaces 
was introduced by Wilson [33]. For this discretization there was also an attempt to prove conver- 
gence of period matrices to their continuous counterparts; see Remark 6.6. Yet another linear the- 
ory for discrete holomorphic functions on regular triangle lattices, due to Dynnikov-Novikov [12], 
was motivated by the theory of integrable systems. A discretization of the Riemann-Roch theorem 
for graphs was introduced by Norine-Becker [1]. 
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Important nonlinear discrete models of complex analysis and Riemann surfaces involve circle 
packings or, more generally, patterns of intersecting or disjoint circles [32, 23, 31]. The linear 
theory of discrete holomorphic functions on quad-graphs can be obtained by linearization from the 
nonlinear circle pattern theory. Discrete holomorphic functions describe infinitesimal deformations 
of circle patterns [4]. A strikingly simple concept of discrete conformal equivalence for triangle 
meshes leads to a new nonlinear model of discrete Riemann surfaces [30], [5]. 

The problem of convergence of discrete analytic functions to their continuous counterparts 
is certainly one of the most important issues. Convergence for the linear theory was proved by 
Courant-Friedrichs-Lewy [8] for square lattices, by Chelkak-Smirnov [7] for rhombic lattices, and 
by one of the authors [25] for orthogonal quadrilateral lattices. Although there are some striking 
convergence results [23, 24, 14], the situation is far from satisfying. 

In this paper we prove convergence of discrete period matrices (Theorem 2.5) and discrete 
Abelian integrals (Theorem 2.6) to their continuous counterparts. In particular, we obtain a 
geometric error estimate (2) for period matrices in terms of circumcircles of the triangulation. We 
also prove a discrete analogue of the Riemann-Roch theorem (Theorem 5.6). 

The proofs use energy estimates inspired by electrical networks, extending those of [25]. En- 
ergy estimates allow to prove convergence of discrete period matrices directly, without proving 
convergence of Abelian integrals first (as in previous approaches). An elementary introduction to 
our method can be found in [26, 27]. 

Main results of the paper are stated in Section 2 and proved in Sections 3-4. Discrete Riemann- 
Roch theorem is stated and proved in Section 5. A generalization of our setup and some open 
problems are given in Section 6. 



2 Main ideas 

2.1 Discrete harmonic and discrete analytic functions 

Let TZ be an oriented 2-dimensional manifold without boundary equipped with a piecewise flat 
metric having isolated conical singularities. The metric of the surface TZ gives rise to a natural 
complex structure (e.g., see [2, Section 1.1.3]). Vice versa, any complex structure on the surface 
can be realized by a piecewise flat metric. In some cases the latter can be easily constructed: for 
instance, for an algebraic curve in 2-dimensional complex projective space the required piecewise 
flat metric has singularities exactly at the branch points of a projection to a 1-dimensional subspace. 

Let T be a geodesic triangulation of TZ such that all the singular points of the metric are vertices 
of T. Denote by T°, T 1 , T 1 , T 2 the sets of vertices, edges, oriented edges, faces, respectively. 
Introduce cotan edge weights by the formula 

, s 1 1 

u{e) = - cot a e + - cot [3 e , 

where a e and (3 e are the angles opposite to an edge e G T 1 in the two triangles sharing the edge e; 
see Figure 1. 

A function u : T° — > M is discrete harmonic, if for each vertex z G T° we have 

^2 v(xy)(u(x)-u(y)) = 0. 

xyf^T 1 : x=z 

Example 2.1. (See Figure 2 to the left.) Let T be a square grid with the SW-NE diagonal drawn 
in each square. Then the sides of the squares have unit weights and the weights of diagonals 
vanish. Thus a function on the vertices of T is discrete harmonic, if the value at a vertex equals 
the average of the values at the neighbors (diagonal neighbors not taken into account). 

For an oriented edge e G T 1 denote by h e G T°, t e G T°, l e G T 2 , r e G T 2 the head, the 
tail, the left shore, the right shore of e, respectively; see Figure 1. Two functions u: T° — > K and 
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Figure 1: Notation associated with an edge e G T 1 ; 



see Section 2.1 for explanation. 



v : T 2 — >■ R are conjugate, if for each oriented edge e G T 



we have 



v(/ e ) - v(r e ) 



u(e)(u(h e ) - u(t e )). 



(1) 



The pair / = (u : T° — > R, v : T 2 — >■ R) of two conjugate functions is called a discrete analytic 
function. We write Re/ := w and Im/ := i>. We write / = const, if both u and v are constant 
functions, not necessarily equal to each other. A direct checking shows that on simply-connected 
surfaces 1Z discrete harmonic functions are precisely the real parts of discrete analytic functions. 

2.2 Discrete Abelian integrals of the 1st kind 

We are going to consider multi-valued functions on T. Informally, a multi-valued function is the 
one whose value changes after performing a nontrivial loop on the surface 1Z. For formal definition 
we need the following notation. 

In what follows assume that the surface 7Z is closed and has genus g > 1 . Denote by p : 7Z — > 7Z 
the universal covering of 1Z and by p: T — > T the induced universal covering of T. Fix closed 
paths «i, . . . , a g , /3\, . . . , (3 g : [0, 1] — > 1Z forming a standard basis of the fundamental group iii(7Z). 
Each closed path a G {a±, . . . ,a g , (3\, . . . , (3 g } determines the deck transformation t a : 1Z — > 1Z, 
i.e., the homeomorphism such that p o t a = p and t a (a(Q)) = 5(1), where 5: E — > 1Z is a lift of 
a : [0, 1] — > 1Z. The induced deck transformation of T is also denoted by t a : T — > T ■ 

A multi-valued function with periods A\, . . . , A g , B%, . . . , B g G C is a pair of functions / = 
(Re/ : T° — > R, Im/ : T 2 — > M) such that for each k = 1, . . . , g and each z G T°, if G T 2 we have 



The numbers Al, . . . , A g and B±, . . . , B g are called the A-periods and the B-periods of the multi- 
valued function /, respectively. A multi- valued discrete analytic function is called a discrete Abelian 
integral of the 1st kind. 

Example 2.2. (See Figure 2.) Let 1Z be the surface obtained by gluing opposite sides of a square. 
Let T be the "triangulation" obtained by drawing the SE-NW diagonal of the square. Then 1Z 
is the plane and T is the infinite square grid in the plane with the SE-NW diagonal drawn in 
each square. The deck transformations t ai ,t^ : 7Z — > 7Z are unit translations along the sides of the 
squares. For a face w G T 2 denote by w* the circumcenter of w. Then a discrete Abelian integral 
of the 1st kind with the period A\ G C is given by the formulae Kef(z) := Ke(Aiz) for each vertex 
z G T° and Imf(w) := lm(Aiw*) for each face w G T 2 . Actually these formulae produce discrete 
Abelian integrals of the 1st kind on any triangulated flat torus 1Z = C/(Z + Z77), where rj G C — R. 



Re/(t Qfc z)-Re/(z) 
Im/(t afc w) - Im/(u>) 



ReA fc ; 



ae/(t /9 ^)-Re/(2) = ReB Jfe ; 
Imf(tp h w) - lmf(w) = lmB k . 
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r A r = s 1 x s 1 

Figure 2: The universal covering of a torus; see Example 2.2 for the details. 

2.3 Discrete period matrices and their convergence 

To proceed we need the following theorem essentially due to Mercat. 

First Existence and Uniqueness Theorem 2.3. For any numbers A±, . . . , A g G C there exist 
a discrete Abelian integral of the 1st kind with A-periods A l7 . . . , A g . It is unique up to constant. 

For each I = l,...,g denote by (f) l T = (Re<f) l T : T° -»■ E,lm0^: T 2 -»■ R) the unique (up to 
constant) discrete Abelian integral of the 1st kind with A-periods given by the formula Ak = 5m, 
where k = 1, . . . ,g. The g x g matrix II7- whose l-th column is formed by the B-periods of <fr l <j-, 
where I = 1, . . . ,g, is called the period matrix of the triangulation T. It is a discrete counterpart 
of the period matrix LT^ of the surface 1Z. For a g x g matrix II denote ||II|| := ySi<j^Kp |IIw| 2 - 

Example 2.4. Let T be an arbitrary triangulation of a flat torus TZ — C/ (Z + Zrj) with rj G C — R. 
Then II7- — rj — Tin by the explicit construction of discrete Abelian integrals in Example 2.2. 

Let us state the main result of the paper. 

Convergence Theorem for Period Matrices 2.5. There are two constants Const^, const^ > 
depending only on the metric of the surface 1Z such that for any triangulation ToflZ with maximal 
circumradius of triangles r < const^ we have 

||n r -n^|| < Const^-r. (2) 

In particular, this implies that for a sequence of triangulations of the surface TZ such that the 
angles of all triangles are bounded from ir and the maximal edge length of n-th triangulation tends 
to zero the discrete period matrices converge to the period matrix LT^. 

Theorems 2.3 and 2.5 are proved in Sections 3.3 and 3.7, respectively. 

2.4 Convergence of discrete Abelian integrals of the first kind 

To state the next result we need the following notions. We say that a discrete Abelian integral 
(f) l T = (Re <fij- : T° — > R, Im <$ r : T 2 — > R) of the first kind is normalized at a vertex z G T° and a 
face w G T 2 , if Re(j) l T (z) = lm.<pj-( w ) — 0. Similarly, denote by 4> l n \ TZ — > C the Abelian integral 
of the first kind with the A-periods given by the formula Ak = 5m, where k = 1, . . . , g. We say 
that (j) l n : TZ — > C is normalized at a point z G TZ, if 4>ti(z) = 0. 

Let {Tn} be a sequence of triangulations of the surface TZ such that the singular points of the 
metric of TZ are vertices of each triangulation T n (the piecewise-flat metric on the surface TZ is fixed 
and does not depend on n). A triangulation T n is Delauney, if for each edge the sum of opposite 
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angles in the two triangles containing the edge is less than tt. The sequence of triangulations {T n } 
is nondegenerate uniform, if there is a constant Const (not depending on n) such that for each 
member of the sequence: 

(D) for each edge the sum of opposite angles in the two triangles containing the edge is less than 
7r — 1 /Const (in particular, the triangulation is Delauney). 

(U) the number of vertices in an arbitrary intrinsic disk of radius equal to the maximal edge 
length is less than Const. 

A sequence of pairs of functions /„ = (Re/ n : 7^° — > R, Im/ n : 7^ 2 — > R) converges to a function 
/ : TZ — > C uniformly on each compact subset, if for each compact set K C TZ we have 

max \Ref n (z) — Kef(z)\ —> and jnax \lmf n (xyz) — Irnf(2:)| — > as n — > oo. 

Convergence Theorem for Abelian Integrals 2.6. Let {T n } be a nondegenerate uniform 
sequence of Delauney triangulations of TZ with maximal edge length approaching zero as n — >■ oo. 
Let z n G 7^° be a sequence of vertices converging to a point z G 1Z and w n G 7^ 2 be a sequence 
of faces such that z n G w n . Then for each 1 < I < g the discrete Abelian integrals (p l Tn = 
(Re (f)j- : — > R, Im cf) l Tn : — > R) o/ i/ie first kind normalized at z n and w n converge to the 
Abelian integral 4> l n : 1Z — > C of the first kind normalized at Zq uniformly on each compact subset. 

This theorem is proved in Section 4.5. 

3 Convergence of discrete period matrices 
3.1 Riemann bilinear identity 

We say that u: T° — > R is a multi-valued function with periods Ax, . . . , A g , Bx, . . . ,B g G R, if 
for each k = 1, ... ,g and each z G T° we have uit ak z) — u{z) = Ak and u{tp k z) — u(z) = B^. 
A multi-valued function v : T 2 — > R is defined analogously. In what follows identify T°,T l ,T 2 
with some subsets of T°,T 1 ,T 2 which project bijectively to T°,T\T 2 , respectively. Introduce a 
bilinear product of two multi- valued functions u : T° — >• R and v : T 2 — > R by the formula 




eeT 1 



where each edge e G T 1 is oriented arbitrarily. 



Riemann Bilinear Identity 3.1. Let u: T° — > R and m': 7^ — > R &e two multi-valued functions 
with periods Ax, . . . , A g , Bx, . . . , B g and A[, . . . , A' , B[, . . . , B' g , respectively. Then 
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(u,u') = J2(A k B' k -B k A' k ). 



k=i 



For the proof we need the following auxiliary assertion. 



Stokes' Formula 3.2. Let T C 7" be a triangulated polygon. Then 



HQ - v{r e )){u{h e ) - u{t e )) = < l e)HQ - u(h e )) 




where each edge e G (dT) 1 is oriented so that l e C T. 
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Proof. In the case when the polygon T is a single triangle of the triangulation T Stokes' formula 
is checked directly. Sum such Stokes' formulas for all triangles of the triangulated polygon T. 
Each boundary edge e contributes the value v(l e )(u(t e ) — u(h e )) to the right-hand side of the sum. 
Each interior edge e appears once with each of the two possible orientations and hence contributes 
(y(l e ) — v(r e ))(u(t e ) — u{h e )) to the right-hand side of the sum. The left-hand side of the sum 
vanishes. Moving the terms (v(l e ) — v{r e )){u{t e ) —u{h e )) to the left-hand side, we get the required 
formula. □ 

Proof of Riemann Bilinear Identity 3.1. Assume without loss of generality that the images of the 
paths ax, . . . ,a g , fti, . . . , (3 g are contained in the union of the edges of T. Cut the surface 1Z along 
these paths. We get a polygon T with 4g sides. Consider T as a subset of the universal covering T. 

Denote by Si, 0i, Sj" 1 , /3-f 1 , . . . , a g , (3 g , S" 1 , the sides of the polygon T in the counterclock- 
wise order, so that they cover the paths «i, (3i, a^ 1 , j3~ , . . . , a g , (3 g , a" 1 , f3~ 1 , respectively. Denote 

~df = Si U 0i U • • ■ U a g U %. 

Identify T 1 with T 1 U (^dT) 1 . Orient all the edges e of dT so that l e C T . 
By Stokes' Formula 3.2 we get 

(it, u) = £ (u'(l e )-u'(r e ))(u(h e )-u(t e )) 

= E u\l e ){u{t e )-u{h e ))+ (u'(l e )-u'(r e ))(u(h e )-u(t e )) 
= E u(l e )(u{t e ) - u(h e )) - E u'{l e ){u{t e ) - u{h e )) 

Since dT is the union of the sides Si, /3i, Sjf 1 , /Jf 1 , . . . a g , /3 g , S" 1 , f3~ l this equals 



= E E AQW**) - U{h e )) - E «U)(u(te) - 
+ E( E w'(^)K*e)-w(^e))- E U\l e ){u{t e )-U{h e )) 

Since u' : T 2 — >• K is multi- valued with periods A[, . . . , A' , B[, . . . , B' g this equals 

9 9 

Since the deck transformations t ak : T — > T (respectively, tp k : T —> T) takes the head of the 
oriented segment S^T 1 (respectively, (5 k l ) to the tail, the sequence of equalities finishes as follows: 



J2(A k B' k - B k A' k 



k=i 

□ 
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3.2 Energy 

The energy of a multi- valued function u : T° — >■ R is the number 

xy&T 1 

Energy Conservation Principle 3.3. Let f = (Re/: T° — >■ R, Imf: T 2 — >• R) 6e a discrete 
Abelian integral of the 1st kind with periods A 1 , . . . , A g , Bx, . . . , B g . Then 

g 

E T (Ref) = -lmJ2A k B k . 

k=l 

Proof. By identity (1) and the Riemann Bilinear Identity 3.1 we have 
E T (Ref) = v(e)(Ref(h e ) - Re/(t e )) 2 

eGT 1 

= J2(Imf(l e ) - lmf(r e ))(Ref(h e ) - Ref(t e )) 

eeT 1 

= (Re/,Irn/} 

9 

= ^2(ReA k lmB k - ReB k lmA k ) 

k=l 

9 

= -1m^2 A k B k , 

k=l 

where each edge e G T 1 is oriented arbitrarily. □ 

Energy Positivity Principle 3.4. For each nonconstant multi-valued function u: T° — > R we 
have E T {u) > 0. 

Proof. Since the metric of the surface 1Z may have singularities only at the vertices of T it follows 
by [11, §4] that £7- (it) = J2 xyz eT 2 Ixyz \^Ixy Z u\ 2 dA, where I xyz u: xyz — > R is a linear extension of 
it to a triangle xyz G T 2 . This proves the principle. □ 

Variational Principle 3.5. A multi-valued discrete harmonic function has minimal energy among 
all the multi-valued functions with the same periods. 

Proof. Consider the finite-dimensional space of multi-valued functions T° — > R with the same 
periods as a given multi- valued discrete harmonic function u: T° — > R. By Energy Positivity 
Principle 3.4 it follows that the energy E-j-(-) is a (nonstrictly) convex functional on this space. 
Since it: T° — > R is discrete harmonic it follows that it is a critical point of the functional Ej-(-). 
Thus it is a local minimum. □ 

3.3 Existence and Uniqueness Theorem 

Corollary 3.6. Let f = (Re/: T° — > R, Im/: T 2 — > R) fre a discrete Abelian integral of the 1st 
kind with periods A\, . . . , A g , B\, . . . , B g . Then the following 3 conditions are equivalent: 

1. f = const; 

2. Ai = ■ ■ ■ = A g = 0; 

3. ReAi = • • • = ReA g = ReSi = • • • = ReB g = 0. 
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Proof. The implications 1=^2 and 1=^3 are obvious. Let us prove that 2=^1; the implication 3=^1 
is proved analogously. By Energy Conservation Principle 3.3 we have 

9 

E T (Ref) = ^2(ReA k lmB k - ReB k lmA k ) = 0. 

k=l 

Hence by Energy Positivity Principle 3.4 we get Re/ = const. By definition of a discrete analytic 
function we get also Im/ = const. Hence / = const. □ 

Proof of First Existence and Uniqueness Theorem 2.3. Uniqueness. Let /, /' be two discrete Abelian 
integrals of the 1st kind as required. Then / — /' is a discrete Abelian integral of the 1st kind with 
vanishing A-periods. By Corollary 3.6(2=^1) it follows that f — f = const. 

Existence. Denote by / = (Re/ : T° — > R, Im/ : T 2 — > R) an unknown multi- valued function 
with real and imaginary parts of the A-periods A\, . . . , A g considered as parameters. 

Consider the following system of linear equations. The variables are the values of the function 
/ = (Re/ : T° -»■ R,Im/: T 2 -»■ R) at the elements of the sets T° C T°, T 2 C T 2 , and also 
unknown real and imaginary parts of the B-periods B\, . . . , B g of this function. Express the 
values of the function / = (Re/ : T° — > R, Im/ : T 2 — >■ R) at the vertices of T° not belonging to 
T° C T° through the values at T° and the periods A±, . . . ,A g , Bi, . . . , B g . Similarly, express the 
values of the function / = (Re/ : T° — > R, Im/ : T 2 — )■ R) at the faces of T 2 not belonging to 
T 2 C T 2 through the values at T 2 and the periods A±, . . . ,A g , B x , . . . , B g . For each edge e G T 1 
choose an orientation and write one linear equation (1). Also pick up a vertex z G T°, a face 
w G T 2 , and write two equations Ref(z) = 0, Im/(w) = 0. 

Denote by \F\ the number of elements in a finite set F. We have written a system of IT 1 ] + 2 
linear equations in |T°| + IT 2 ! + 2g variables. By the Euler formula the number of equations in the 
system equals the number of variables. The free terms of these equations are linear combinations 
of the parameters ReAi, ImAi, . . . , KeA g , lmA g . By uniqueness part of the theorem it follows that 
the homogeneous system corresponding to A\ — ■ ■ ■ = A g = has only trivial solution. Then 
by the finite-dimensional Fredholm alternative it follows that for arbitrary A\,...,A g G C the 
nonhomogeneous system has at least one solution. Thus there exists a discrete Abelian integral of 
the 1st kind as required. □ 

The following theorem is proved analogously. 

Second Existence & Uniqueness Theorem 3.7. For each P = (Ax, . . . , A g , Bi, . . . , B g ) G R 2g 

there is a unique (up to constant) discrete Abelian integral 4>t,p = (Re0r,P : 7~° — > R, lm<pj- t p: T 2 — > 
R) of the 1st kind whose periods have real parts A±, . . . , A g , B\, . . . , B g , respectively. 

3.4 Discrete period matrices 

For each I = 1, . . . ,g denote by (j) l T , = (Re0^* : T° —> R, Im^, : T 2 — > R) the unique (up to 
constant) Abelian integral of the 1st kind with A-periods given by the formula A k = i5 k i, where 
k = 1, . . . ,g. The g x g matrix n-p. whose l-th column are the B-periods of <p l T , divided by % is 
called the dual period matrix of the triangulation T . 

Example 3.8. The map taking the vector of A-periods of a discrete Abelian integral of the first 
kind to the vector of its B-periods is not necessarily C-linear, and thus n^* ^ H^ in general. For 
instance, if T is obtained from the side surface of a regular square pyramid by identifying the 
opposite sides of the base then Uj- = 2i/y/3, U-r* = i\/3/2. 

Period Matrix Lemma 3.9. The matrices Imllx and Imrif. are symmetric positively definite, 
and Ren r » = (Ren r ) T . 
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Remark 3.10. The discrete period matrix by Mercat [18] is IIg := (n7- + n-p.)/2. The matrix Ilg is 
symmetric with positively definite imaginary part. Our proof of Convergence Theorem for Period 
Matrices implies an analogue of error estimate (2) for this matrix as well: ||IIg — II^H < Const-^-r. 

For the proof of the lemma we need the following auxiliary assertion. 

Lemma 3.11. Let f = (Re/: f° -»■ R,Im/: f 2 R) and f = (Re/': f° -»■ R,Imf : f 2 -»■ R) 
6e two discrete Abelian integrals of the first kind with periods Ax, . . . , A g , B\, . . . , B g and A[, . . . , A' g , 
B[, . . . , B' g , respectively. Then 

Im j^{A k B' k - B k A' k ) = 0. 
fc=i 

Proof. By Riemann Bilinear Identity 3.1 we have 

9 9 

Im ^2(A k B' k - B k A' k ) = ^2(ReA k lmB' k - ReB k lmA' k - ReA' k lmB k + ReB' k lmA k ) 

k=l k=l 

= (Re/,Im/')-(Re/',Im/) 

= J2 (W(/ e ) - Im/'(r e ))(Re/(/ ie ) - Re/(t e )) 

eGT 1 

- (W(4) - Im/(r e ))(Re/ / (/i e ) - Ref'(t e )) 

eer 1 
= 0, 

where each edge e G T 1 is oriented arbitrarily, because by formula (1) the respective summands 
in the latter two sums cancel each other. □ 

Proof of Period Matrix Lemma 3.9. Applying Lemma 3.11 for / := f := (fxL we get (Imll7-)y— 
(Imllj-)ji = 0, that is, Imll7- is symmetric. Applying Lemma 3.11 for / := (f)j-, f '■= <$y* we get 
(Ren r *) - - (Ren r )j< = that is, R e n r * = (R e n r ) T . 

Finally, let / be a discrete Abelian integral of the 1st kind whose A-periods A\, . . . , A g are 
arbitrary real numbers not vanishing simultaneously. Then the B-periods of this integral are 
equal to Bi = Ylk=i(^T)kiA k - Hence by Energy Conservation Principle 3.3 and Energy Positivity 
Principle 3.4 it follows that 

9 

£ (hnn. T )kiA k Ai = - M^kB k ) = E r (Ref) > 0. 

l<k,l<g k=l 

Thus IH1II7- is positively definite. Analogously, IH1II7-. is symmetric and positively definite. □ 

Denote ur,p := Re^r.P; where <Pt,p is th e discrete Abelian integral of the 1st kind defined in 
Second Existence and Uniqueness Theorem 3.7 for each vector P G R 29 . 

Lemma 3.12. The energy E-j-(ut,p) is a quadratic form in the vector P £ R 29 with the block 
matrix 

_/Ren r »(Imn r -1 Ren r + Imn r (ImrirO^RelLA 
hr '-~y Renr^Imllr*)- 1 (Imn r ,) _1 )' 

Proof. Denote by Ax, ■ ■ ■ , A g , B\ . . . , B g £ C the periods of the discrete Abelian integral 4>t,p 
of the 1st kind. Denote A := (Ax,...,A g ), B := (Bx,...,B g ). Then P = (Re A, ReB) = 
(ReAx, • • • , ReA g , ReBi . . . , ReB g ). By First Existence and Uniqueness Theorem 2.3 we get 

cf) T P = ReAx </>- H h ReA g cj) 9 T + ImAx (j) l T * H h lmA g (j) 9 T * + const. 
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Thus 

ReB = ReU T ReA - Imn T *ImA; 
ImB = Imn r ReA + ReLVJmA 

Solving in ImA and ImB we get 

1mA = (Imn r ») _1 Ren r ReA - (ImLV^Re-B; 

ImB = (lmU T + Renr^Imnr*)" 1 Ren r ) ReA - Ren r »(Imn r »)- 1 ReS. 

Applying Energy Conservation Principle 3.3 for / := <Pt,p we get the required expression for the 
matrix of Ej-(ut^p). □ 



3.5 Energy in continuous theory 

We need the following result from continuous theory proved analogously to Lemma 3.12. Let 
P = (A\, . . . , A g , Bi, . . . , B g ) G M 29 and let 4>n.p : 7?. — ?■ C be an Abelian integral of the first kind 
whose periods have real parts A\, . . . , A g , B\, . . . , B g . Denote u-n,p '■= Re07e,p. The (Dirichlet) 
energy of a multi- valued continuous piecewise-smooth function u : 7Z — > R is 

E n {u) := / \Vu\ 2 dA. 
Jn 

(Here the gradient Vw is well-defined on 1Z and not just on 1Z.) 

Lemma 3.13. The energy Eh(uti,p) is a quadratic form in the vector P G M? 9 with the block 
matrix 

/Ren^(imn^)- 1 Ren^ + imn 7e (imn^)- 1 R e n^\ 
71 ■ V Ren^(imn^)- 1 (imn K )^ )■ 

3.6 Convergence of energy 

In this subsection we prove the following result. 

Energy Convergence Lemma 3.14. For any vector P G M 29 there exist two constants Const p^, 
const > depending only on P and 1Z such that for any triangulation T of 1Z with maximal 
circumradius of triangles r < const p.n we have 

\Et{ut,p) ~ E n (u n ,p)\ < Const • r. 

For a function u: T° — > R define its interpolation to be the linear extension lq-u: 7Z — > R to 
each triangle of T. 

Lemma 3.15. For any function u: T~° — > R we have Eq-(u) = E-jiilq-u). 

Proof. Since the metric of the surface 1Z may have singularities only at the vertices of T it follows 
by [, 1, §4] that Ef{u) = J^weT 2 I \^lTu\ 2 dA, and the lemma follows. □ 

Lemma 3.16. For any P G M 29 we have E^(ur,p) > E^u^p). 

Proof. This follows from Lemma 3.15 and the continuous variational principle because the inter- 
polation /7-M7-P : TZ — > R is continuous. □ 

The following lemmas are analogous to [25, Claims 4.2, 4.5, and Lemma 2.3]. 
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Projection Lemma 3.17. For any triangle xyz C R 2 and any vector it el 2 we have 

v < 2 secfmaxjZx?/;?, Zyzx, Zzxy}) ■ — j — j 1 j — j 1 j — j — . 

\ \xy\ \yz\ \zx\ J 

Proof. Denote a = m&x{Zxyz, Zyzx, Zzxy}. Then the value /-(xy, zx) is between a and n — a. 
Since Z(xtj, z&) = ±Z(lf, xtj) zx) it follows that at least one of the angles in the right-hand 

side, say, the first one, does not belong to the interval (§ — f > § + f )• Then 

v < sec — ■ — ; — ; — < 2 sec a ■ 



\xy\ \xy\ 



□ 



Gradient Approximation Lemma 3.18. Let xyz C R be a triangle of circumradius r. Let 
u: xyz — > R be a C 2 function. Denote by I X yzU: xyz — > R its linear interpolation in xyz, i.e., the 
linear function such that I xyz u(x) = u(x), I X yzu(y) = u{y), I xyz u(z) = u(z). Then 

max \Vu(w) — VI xyz u(w) \ < 12r max \D 2 u[w)\. 

wdxyz wdxyz 

Proof. By the Rolle theorem there is a point w G xy such that (V«(to) — V I xyz u(w)) -xtj/\xtj\ = 0. 
Since VI xyz u is a constant vector it follows that (Vh — \/I xyz u)-xy/\xy\ < Di&m(xyz) max xyz \D 2 u\ 
in the triangle xyz. The same inequality holds with xy replaced by yz or zx. Thus the lemma 
follows from Projection Lemma 3.17 because Diam(xyz) sec(max{Zxyz, Zyzx, Zzxy}) = 2r. □ 

Lemma 3.19. For any C 2 multi-valued function u: 1Z — > R there exist two constants Const Uj 7^. ; 
const,^ > depending only on u and 1Z such that for any triangulation T of 1Z with maximal 
circumradius of triangles r < const u n we have 



\E T (u \ f ) - E n (u)\ < Const^ 



r. 



Proof. Since 1Z is compact and |Vw|, \D 2 u \ : 1Z — > R are continuous functions it follows that \Vu\ 
and \D 2 u\ are bounded. Then by Lemma 3.16 and Gradient Approximation Lemma 3.18 we get 

\E T {u \ f0 ) - E n {u)\ = \E n {I T u \ f0 ) - E n {u)\ 

^2 / (2Vu + VI xyz u- Vti)(Vtt- VI xyz u)dA\ 

r-j — 9 J xyz 



xyzdT 2 



< (2 max I Vm| + 12r max \D 2 u\) ■ 12r max \D 2 u\ ■ Areaft. 
n n n 

The latter expression is a quadratic function q(r) in r. It remains to set const u ^ := 1, Const^ := 
?(!)■ ^ □ 

Proof of Energy Convergence Lemma 3.14- The Lemma follows from the following sequence of es- 
timates provided by Lemma 3.16, Variational Principle 3.5, and Lemma 3.19: 

< E T (u TjP ) - E n {u nt p) < E T (u n ^ P |^ ) - E n (u n ^ P ) < Const • r. 

□ 



3.7 Convergence of discrete period matrices 

Proof of Convergence Theorem for Period Matrices 2.5. By Lemmas 3.12 and 3.13 the energies 
Eq-(uq-^p) and E^u^p) are quadratic forms in P G R 29 with the block matrices Ej- and En, 
respectively. Thus by Lemma 3.14 there exist two constants Const-^, const-^ > depending only 
on 1Z such that for any triangulation T of 1Z with maximal circumradius of triangles r < const^ 
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we have \\Eq — En\\ < Const-^ • r. In particular, ||(ImIl7-») _1 — (Im n^) — 1 1| < Const^ • r for each 
r < const-fc. Hence there exist two new constants Const^, const^ > depending only on 71 such 
that || (Im II7-*) || < Const^ for each r < const^. Assume further that r < const^ and r < const^. 
By the inequality \\Ej — En\\ < Const^ • r and the properties of the norm we have 



Const^-r > ||(Imn r ») _1 ReII r - (Imn^)" 1 R e n^|| 

= iKimnrO-^Rerir - R e n^) - ((imn^)- 1 - (imnr-)- 1 ) Ren^|| 

> (KimnrOir 1 • ||Ren r -Ren^|| - iKimn^^-^mnr.) -1 !! • ||Ren*| 

> (Const^) -1 ■ II RellT — RelT^H — Const-^. • r • || ReH^H 



Thus || RelT? — ReH^H < Const^ • r, where Const^ := Const^ • Const-^. • (1 + || Reil^H). 
Finally, since \\Eq — En\\ < Const^ • r it follows that 

|| Re II r * (Im II r * y 1 Re Tl T + Im Tl T - Re IL^Im n^)" 1 Re - Im || < Const^ ■ r. 

Using the result of the previous paragraph and similar estimates we conclude that there exist a 
new constant Const^ such that || Ren7-*(lmll7-»)~ 1 ReIT7- — Ren-T^ImlX^) -1 ReH^H < Const^-r. 
Thus || IhiiIt — Imil^H < (Const^ + Const^) • r, which completes the proof of the theorem. □ 

4 Convergence of discrete Abelian integrals 

In this section we prove Convergence Theorem for Abelian Integrals 2.6. The proof goes along the 
same lines as the proof of [25, Convergence Theorem 1.2]. It uses several lemmas stated in the 
next few subsections. 

4.1 Convergence of energy 

Lemma 4.1. Let {T n } be a nondegenerate uniform sequence of Delauney triangulations of 1Z with 
maximal edge length approaching zero as n — > 00. Let P n G M 29 be a sequence of 2g- dimensional 
real vectors converging to a vector P G M? 9 . Then Ej- n (uj- n! p n ) —> Etz{uk,p) as n — )■ 00. 

Proof. Both Ej- n (uj- n ,p) and En(un^p) are quadratic forms in P G M? 9 . By Energy Convergence 
Lemma 3.14 the former sequence of forms converges to the latter form. This implies the lemma. □ 

Assume that v[e) 7^ for each edge e G T 1 . Then the energy of a multi- valued function 
v: T 2 — > K. is E-y{v) := J2 e€T i(v(l e ) — t>(r e )) 2 /V(e), where each edge e G T 1 is oriented arbitrarily. 

Lemma 4.2. Let f = (Re/: T° — > R, Im/: T 2 — > R) be a discrete Abelian integral of the 1st kind. 
Then E r (Ref) = E T (lmf). 

Proof. This follows immediately from identity (1). □ 

4.2 Equicontinuity 

The eccentricity of a Delauney triangulation T is the infimum of the numbers Const such that the 
triangulation satisfies conditions (D) and (U) from Section 2.4. 

Equicontinuity Lemma 4.3. Let T be a Delauney triangulation of a planar polygon T , K be 
a compact set inside &T , and u: T° — > R be a discrete harmonic function. Denote by r : = 
min zg ^ Dist(z, &T) and by e the eccentricity of T ■ Then there is a constant Const x, r ,e depending 
only on K, r, e (but not on T, u, z, w) such that 

\u{z) - u{w)\ < Const ^ r , e • E T {u) l/2 ■ ln~ 1/2 (\ + W ^ J ( 3 ) 

for any z,w G K fl T° . 
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Remark 4.4. This result remains true, if one replaces a planar polygon T by a piece of a cone with 
the aperture 9, and the constant Const^e by a constant Const^^e depending also on 6. The 
generalization is reduced to the previous result by projection of the cone into the plane along the 
axis. 

Proof. The lemma is proved analogously to a similar estimate for quadrilateral lattices [25, Equicon- 
tinuity Lemma 2.4]. Actually one can derive the lemma from that estimate as follows; cf. [25, 
Section 5.1]. In each triangle of T, draw 3 segments joining the circumcenter with the vertices. 
Erase hanging edges from the obtained graph. We get an orthogonal quadrilateral lattice Q; see 
Figure 1. The function u: T° — > R is the restriction of a discrete harmonic function Q° — > R with 
the same energy. By [25, Equicontinuity Lemma 2.4] the lemma follows. □ 

4.3 Harmonicity of a uniform limit 

Let C C be a planar domain. A sequence of triangulated polygons {T n } approximates the domain 
Q, if for n — > oo: 

• the maximal distance from a point of dT n to the curve dVt tends to zero; 

• the maximal edge length of the triangulation T n tends to zero. 

A sequence of functions u n : T® — > R converges to a function u : Q — > R uniformly on each compact 
subset, if for each compact set K C Q we have 

max \u n (z) — u(z)\ — > as n — > oo. 

Lemma 4.5. Let {T n } be a nondegenerate uniform sequence of Delauney triangulations of polygons 
T n approximating a bounded domain flcC. Let u n '- T® — >■ R be a sequence of discrete harmonic 
functions converging to a continuous function u\ Q — > R uniformly on each compact subset. Then 
the function u : Q — > R is harmonic. 

This lemma immediately reduces to its counterpart for quadrilateral lattices [25, Lemma 4.15] 
analogously to the proof of Equicontinuity Lemma 4.3 above. 

4.4 Convergence of the real parts 

We say that the function uq- t p : T° — > R (respectively, w^p : 1Z — > R) is normalized at a point z if 
u t,p( z ) — (respectively, un,p{z) = 0). 

Theorem 4.6. Let {T n } be a nondegenerate uniform sequence of Delauney triangulations of 1Z with 
maximal edge length approaching zero asn->oo. Let z n £ T® be a sequence of vertices converging 
to a point z G 1Z. Let P n G R 2fl be a sequence of 2g -dimensional real vectors converging to a 
vector P G R 29 . Then the function uj- nt p n : 7^° — > R normalized at z n converges to u-jip: TZ — > R 
normalized at zq uniformly on each compact subset. 

Proof of Theorem 4-6. Take an arbitrary subsequence {T nk } of the given sequence {T n }. For 
brevity denote 7X := T n . k - Take a sequence of compact sets K\ C K 2 C • • • C TZ such that 
TZ = U'jL^j. Assume that K 1 contains all the points of the converging sequence {zk}. 

We start with a "handle decomposition" of the triangulated surface TZ and estimates for each 
individual handle. For each edge e 6 T/ denote by Q e C TZ the union of the edge e (without 
the endpoints) and the interiors of the two triangles of 7i containing the edge e. For each vertex 
z G Ti denote by Q z C TZ the union of the vertex z, all the edges of T\ containing the vertex z 
(without endpoints), and the interiors of all the triangles of T\ containing the vertex z. Then each 
set fl e is isometric to a planar domain and each set Q z is isometric to a piece of a cone. 



13 



Take an edge e G Ty and an arbitrary compact subset K e C Vt e . Let Tk be the union of all the 
triangles of the triangulation Tk contained in the domain Q e . Let us estimate the right-hand side of 
inequality (3) from Equicontinuity Lemma 4.3 for K := K e and u := wr fc ,P fc fo ■ Since the maximal 

edge length of Tk tends to zero it follows that there is a number k\ such that for each k > k\ the set 
K e is contained inside dTk and Dist(K e , dTk) is bounded from zero. By Lemma 4.1 the sequence 
of energies Eq- k (uj- k ^p k ) is bounded. Thus by Equicontinuity Lemma 4.3 it follows that the function 

u T k ,p k K e nf° nas un tf orm ^y hounded differences, i.e., there is a constant Const not depending on 

k, z, w such that for each k > k\ and z, w G K e n T k ° we have \u% t p h (z) — Uf k ,p k {w)\ < Const. By 

the same lemma the sequence {u-y k ,p k K e nf° } * s equicontinuous, i. e., there is a positive function 

5(e) not depending on k, z, w such that for each k > k\ and z, w G K e D T k ° with \z — w\ < 8(e) 
we have \u Tk)Pk (z) - u TkA (w)J < e. 

Now take a vertex z G 7^° and a compact subset K z C f2 2 . Analogously to the previous 
paragraph using Remark 4.4 we obtain that {ur k ,p k K z nf° } ^ as un if° rm ly bounded differences and 
is equicontinuous. 

Let us combine the obtained estimates. Since the compact set K\ is contained in the union 
of finitely many sets of the form K e and K z it follows that the sequence {ur k ,p k Kxnf } a ^ so 
has uniformly bounded differences and is equicontinuous. Moreover, the sequence is uniformly 
bounded, because all z k G K\ and uj- kj p k (z k ) = 0. Then by the Arzela-Ascoli theorem it follows 
that there is a continuous function U\ : Ky — > R and a subsequence {Ik} of the sequence 1,2,... 
such that li = k\ and uj uPl converges to U\ uniformly in K\. 

Proceed to the next compact set K2. Analogously, there is a function u<i~- K2 — > K and a 
subsequence {mk} of the sequence {Ik} such that mi = l±, m 2 = / 2 , and ur m ,p m converges to w 2 
uniformly on K 2 . Clearly, U\ = u 2 on K\. Thus the extension can be continued, and eventually 
we get a continuous function u: 1Z — > R and a subsequence {p k } of the sequence 1, 2, . . . such that 
u-j- p converges to u uniformly on each compact subset of 1Z. 

Clearly, u: TZ — > R has the same periods P as un,p and u(z ) = 0. By Lemma 4.5 applied 
to each domain Q e it follows that the limit function u: TZ — > R is harmonic in 7?. possibly except 
the vertices of the triangulation Ti (where the metric of TZ may have conical singularities). By 
singularity removal theorem it follows that the continuous function u : TZ — > R is harmonic in the 
whole surface TZ. Thus u = uk,p normalized at z . 

Since the limit function u is unique, it follows that the initial sequence u-j- nt p n : T® — > R, not 
just the subsequence ur k ,p k , converges to un,p uniformly on each compact subset. □ 



4.5 Convergence of discrete Abelian integrals of the first kind 

Proof of Convergence Theorem for Abelian Integrals 2.6. Let P n ,P G M? 9 be the periods of the 
real parts Re (f) l Tn : 7^° — > R and Re <j) l n : TZ — > R, respectively. Then by Second Existence and 
Uniqueness Theorem 3.7 it follows that Re<j) l rn = ur n ,p n and Re0^ = Un,p- By Convergence 
Theorem for Period Matrices 2.5 we have P n — > P as n — >■ 00. Thus by Theorem 4.6 the real parts 
Re(p l Tn : 7^° — > R converge to Re0^: 7?. — >• R uniformly on each compact subset. Convergence of 
the imaginary parts is proved analogously using Lemma 4.2. □ 



5 Discrete Riemann— Roch theorem 

5.1 Discrete Abelian integrals of the second kind 

A discrete Abelian integral of the second kind is an arbitrary multi- valued function / = (Re/: T° — > 
R, Imf : T 2 — > R). If it does not satisfy equation (1) for an edge e G T 1 then we say that / has a 
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simple pole at the edge pe. The value 

res e / := Im/(r e ) - lmf(l e ) + u(e)Ref(t e ) - v(e)Ref(h e ) 

is called the residue of /. In the case when all the periods vanish then the discrete Abelian 
integral of the second kind is called a discrete meromorhic function; it can be viewed as a pair of 
single- valued functions (Re/: T° -> R, Im/: T 2 — > R). 

Third Existence and Uniqueness Theorem 5.1. For each edge e G T 1 there is a unique (up 
to constant) discrete Abelian integral e = (Re0 e : T° — > R, Im0 e : T 2 — > R) of the second kind 
with vanishing A-periods, with the only pole at the edge e, and with res e <p e = 1. 

This theorem is proved analogously to First Existence and Uniqueness Theorem 2.3. The only 
difference is that the value res e e should be considered as a parameter instead of the values of 
A-periods. 

Lemma 5.2. The B-periods of <j) e = (Re0 e : T° — > R, Im0 e : T 2 — > R) are 

B x = Re4> l T »(t e ) — Re(f) l T *(h e ) + iRe(f) l T (t e ) — iRe 4> l T (h e ) , where I = 1, . . . ,g. 

Proof. Let A ly . . . , A g , Bi, . . . , B g be the periods of <fi e . Denote by A[, . . . , A' g , B[, . . . , B' g the peri- 
ods of (j) l T . All the A-periods Ai, . . . , A g , A[, . . . , A' g vanish except for A\ = 1. By Riemann Bilinear 
Identity 3.1 we have 

g 

ImBi = ^2(ReA' k lmB k - ReB' k lmA k - ReA k lmB' k + ReB k lmA' k ) 

k=l 

= (Re0^-,lm0 e ) - (Re(p e ,lm^ T ) 

= (Im0 e (//) - Im0 e (r / ))(Re0J r (^/) - 
/er 1 

- ^ (Im^(Z/) - Im^(r / ))(Re0 e (/ i/ ) - Re0 e (^)) 
/er 1 

= ^ re S/ e - (Re0V(*/) ~^M h f)) 
feT 1 

= Re<f) l T {t e )-Re<f) l T {h e ), 
where each edge / G T 1 is oriented arbitrarily. Analogously Re£>z = Re0 r *(t e ) — Re0^-»(/i e ). □ 

5.2 Discrete Abelian differentials 

The differential of a discrete Abelian integral / = (Re/ : T° — > R, Im/ : T 2 — )■ R) of the 1st kind 
is the function df : T 1 — > R given by the formula 

d/(e) := Re/(/i e ) - Ref(t e ) (4) 

for each edge e G T 1 . The functions w : T 1 — > R which are differentials of discrete Abelian integrals 
of the 1st kind are called discrete Abelian differentials of the 1st kind. 

Discrete Abelian differentials of the 2nd kind are undefined in our setup. 

A discrete Abelian differential of the third kind is an arbitrary map u: T 1 — > R satisfying 
the condition oj{— e) = —u(e). Its integral J u along an oriented path 7 in the graph T 1 is the 
sum of values on the oriented edges of the path 7. The residue res w u of u at a face w G T 2 is 
the counterclockwise integral along the boundary of the face. Given a collection of closed paths 
ai, . . . , a g , (3i, . . . , (3 g : [0; 1] — > T 1 in the graph T 1 forming a standard basis of 7Ti (72.) , define the 
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real periods of 00 as the integrals along these paths. The real periods of u may depend on the 
choice of the paths. 

Similarly, a discrete Abelian differential u : T 1 — > M. of the third kind can be integrated along 
any chain of oriented triangles such that the consecutive triangles in the chain share a common 
edge. For example, if the chain consists of 2 oriented triangles xyz and yzw then the integral of u 
is by definition u(yz)u(yz). The details of the definition for a general chain are left to the reader. 
The residue of u at a vertex z G T° is the number 

res z u := i v(e)u(e). 

eef 1 :t e =z 

The imaginary periods of u are defined analogously to the real ones as the integrals along certain 
fixed chains of triangles. 

A direct checking shows that u>: T 1 — > K is a discrete Abelian differential of the 1st kind if 
and only if its residues at all faces and all vertices vanish. It is also easy to see that for a discrete 
Abelian differential of the 3rd kind the sum of the residues at all the vertices and all the faces 
vanishes. 

Fourth Existence and Uniqueness Theorem 5.3. For any two distinct vertices z,w G T° 
there is a unique discrete Abelian differential d<fr z , w : T 1 — > M of the third kind with res z d<j) ZjW = 
—res w d(f)z t w = i and with all the other residues and the A-periods vanishing. 

Proof. Uniqueness. Let u, u' be two discrete Abelian differentials of the third kind as required. 
Then all residues of u — u' vanish. Hence u — u' is a discrete Abelian differential of the first kind. 
Since all A-periods of u — u' vanish by Corollary 3.6(2)=>(1) it follows that oj = u' . 

Existence. Let u>: T 1 — > K be an unknown function satisfying the condition u(— e) = —u(e). 
Consider its values at all the edges as variables, and the residue ir at the vertex z as a parameter. 
For each vertex y G T° distinct from w and z write one linear equation res^o; = 0. Write the linear 
equation res 2 o; = ir. (The equality res„,u; = — ir is obtained by summation of the written ones 
and therefore is not written.) For each face x G T 2 but one write one linear equation res^o; = 0. 
(Similar equation for the remaining face is obtained by summation of the written ones.) For each 
k — 1, . . . , g write one linear equation u = meaning that the real A-periods of u vanish. Write 
g analogous linear equations meaning that the imaginary A-periods of u vanish. We get a system 
of |T°| + \ T 2 \ + 2g — 2 equations in \T 1 \ variables. By the Euler formula it follows that the number 
of equations equals the number variables. By the uniqueness part of the theorem it follows that 
for r = the system has only trivial solution. Thus by the finite-dimensional Fredholm alternative 
it follows that for r = 1 the system has a solution. This solution is the required discrete Abelian 
differential of the third kind. □ 

Lemma 5.4. For any discrete meromorphic function f = (Re/: T° — > M, Imf : T 2 — > IR) and any 

two distinct vertices z,w G T° we have 

Ref(z) - Ref(w) = d<j) Z)W (e) ■ res e f 

eef 1 

Proof. The lemma follows from the sequence of equalities 

^2 d^zA 6 ) ■ Ies ef = 5^ d ^^M ' ( lm f( r e) ~ W(U + v(e)Ref(t e ) - u(e)Ref(h e )) 

eef 1 eef 1 

= lm f( x ) ^2 ^.«>( e ) + Re ^) v ( e ) d( t>z,w{e) 
zeT 2 eeft-.r^x y£T° ee fi. te=y 

= ^ Imf(x)res x d(fr ZjW - i ^ Ref(y)res y d(j) ZjW 

x&T 2 y eT° 
= Ref(z)-Ref(w). 

□ 

Similar results hold, if one takes two distinct faces z,w G T 2 instead of two distinct vertices. 
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5.3 Discrete Riemann— Roch theorem 



A divisor on T is an arbitrary map T° U T 1 U T 2 — > {0, ±1}. The divisor (/) : 7~° U T 1 U T 2 — > 
{0, ±1} of a discrete meromorhic function / = (Re/: T° — >■ M, Im/: T 2 — » M) is defined by the 
formula: 

mn f +1 > ifRe/W = 0; f nfe] = [-l, ifres e /^0; mM = f+1, if W (>) = 0; 

\0, ifRe/(*)^0; \o, if resj = 0; J " \o, if Imf(w) / 0; 

for each 2 G T°, e G T 1 , u> G T 2 . The divisor (cj) : T° U T 1 U T 2 — >■ {0, ±1} of a discrete Abelian 
differential u; : T 1 — > M of the third kind is defined by the formula: 

. -1, ifres 2 u;^0; f+1, ifw(e) = 0; f-1, if res^u; ^ 0; 

M( z )= n n M(e) = < n , . (w)(w) = -j 

(J, if res z u; = 0; 10, if ui[e) 0; 10, 11 res^u; = I). 

Remark 5.5. As opposed to the continuous theory, on a discrete Riemann surface divisors do not 
form an Abelian group. This reflects the fact that the product of discrete analytic functions is 
not necessarily a discrete analytic function again. Also the notions of multiple zeroes and multiple 
poles are undefined in our setup. 

The degree deg D of a divisor D : T° U T 1 U T 2 -> {0, ±1} is the sum of its values. The divisor 
L> : T° U T 1 U T 2 -)■ {0, ±1} is greater or egwa/ to a divisor £>' : T° U T 1 U T 2 -)> {0, ±1}, if for each 
z G T° U T 1 U T 2 we have -D(z) > D'(z). Denote by 1(D) the dimension of the space of discrete 
meromorphic functions on T with divisors greater or equal to the divisor D. Denote by i(D) the 
dimension of the space of Abelian differentials of the third kind on T with divisors greater or equal 
to the divisor D. A divisor D is admissible, if D(z) < for each vertex z G T°, D(e) > for each 
edge e G T 1 , and < for each face w G T 2 . 

Riemann— Roch Theorem 5.6. For any admissible divisor D on a triangulated surface of genus 
g we have 

l(—D) = deg-D — 2g + 2 + i(D). (5) 

Remark 5.7. This theorem is a real version of the classical Riemann-Roch theorem: D<r) = 
deg D-ji — g + 1 + i(D n ), where is a divisor on a Riemann surface 1Z and -D^), ^(D-r) are 
the complex dimensions of continuous counterparts of the spaces from the definition of l(—D),i(D). 
The discrete counterpart D of a divisor .D-r. has twice larger degree. An informal explanation is 
that the counterpart of one complex equation describing vanishing of a meromorphic function at 
a point are two real equations describing vanishing of the real and imaginary parts of a discrete 
meromorphic function. 

Proof of Riemann-Roch Theorem 5.6. Let D : T° U T 1 U T 2 — > {0; ±1} be an admissible divisor. 
Denote D^ = D\fi and D = D | jo u js . By Corollary 3.6 a discrete Abelian integral / = 
(Re/ : T° — > ffi, Im/ : T 2 — )■ R) of the 2nd kind with vanishing A-periods whose poles belong to 
the support of the divisor D^ can be written in a unique way as 

/ = 2J ^ e + const 

eGsuppDoo 

with A e = res e / G R. By Lemma 5.2 the B-periods of the discrete Abelian integral / = (Re/ : T° — > 
M, Im/ : T 2 — 7- K) of the 2nd kind vanish if and only if for each k — 1, . . . , g we have 

£ A e #*-(e) = 0; (6) 

eesuppDoo 

A e #^(e) = 0. (7) 

eGsuppDoo 
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Assume that equations (6)— (7) hold. Then by Lemma 5.4 the restriction of / = (Re/: T° — > 
R, Im/ : T 2 — > R) to the set suppDo is a constant if and only if for each pair of elements z,w G 
T° fl suppD or z, w G T 2 fl supp-Do we have 

Kd<j> z , w {e) = 0. (8) 

eSsuppDoo 

In what follows assume that equations (8) are written only for those pairs (z, w) which form certain 
maximal trees in the two complete graphs on the sets T° fl supp_D and T 2 fl supp-D , respectively. 
(The equations for the remaining pairs of vertices follow from the written ones.) 

Consider (6)-(8) as a system of linear equations in variables A e for all e G supply. By the 
previous paragraph the number l(—D) equals to the dimension of the space of solutions of this 
system. Thus the rank of the system (6)-(8) equals deg-Doo — l(—D). 

Transposing the matrix of the system (6)-(8), we get a new system of linear equations. Clearly, 
i(D) equals to the dimension of the space of solutions of the new system. Thus the rank of the 
new system equals deg D + 2g — 2 — i(D). Since transposition does not change the rank we obtain 
the required identity deg-Doo — l(—D) = degD + 2g — 2 — i(D). □ 

Example 5.8. The Riemann-Roch theorem in the particular case of a triangulated flat torus 
implies the following: 

• there are no nonconstant discrete meromorphic functions with exactly one simple pole; 

• there exists a nonconstant discrete meromorphic function with exactly two simple poles at 
two given edges if and only if the edges are parallel to each other. 

This is proved by a direct computation of the right-hand side of (5) using explicit construction of 
discrete Abelian integrals in Example 2.2. 



6 Generalizations and open problems 
6.1 Quad- surfaces 

Let us introduce a generalization of discrete Riemann surfaces considered in this paper above. 
This generalization has a physical interpretation in terms of alternating-current networks; cf. [25, 
Section 5.2]. 

Let Q be a cell decomposition of the surface 71 with quadrilateral faces. Assume that each 
vertex of Q is paint either black or white so that the endpoints of each edge have different colors. 
Assume also that any two faces have at most one common edge. A chart on Q is an orientation- 
preserving isometry of a face and a (not necessarily convex) quadrilateral z\Z%z%z± in the complex 
plane C. An atlas on Q is a collection of charts, one for each face (no agreement of charts for 
distinct faces is assumed). A cell decomposition Q with a fixed atlas is called a quad-surface. 

A function /: Q° — > C is discrete analytic on Q, if the difference quotients along the two 
diagonals of each face are equal, i. e., 

/(fi) - f(z 3 ) f(z 2 ) - /(z 4 ) 

Z\ — Z 3 Z2 — Z4 

for each face z\z 2 z 3 z± G Q 2 [19]. 

A multi-valued function on Q with periods A 1; . . . , A g , B 1; . . . , B 9 , A 1; . . . , A g , Bi, . . . ,M g G C 
is a function / : Q° — > C such that for each k — 1, . . . , g, each black vertex z E Q°, and each white 
vertex w G Q° we have 

f{t ak z) - f(z) = A k ; f(t Pk z) - f{z) = B fc ; 

f(t ak w) - f(w) = A k] f(t h w) - f(w) = B fc . 
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The numbers Ai, . . . , A g , A±, . . . , A g are called the A-periods of the multi-valued function /. A 
multi-valued discrete analytic function is called a discrete Abelian integral of the 1st kind on Q. 
Introduce a skew-symmetric bilinear product of two functions /, /' : Q° — > C by the formula 

(/, /') := £ «M) - /(*»))(/'(**) - /'&)) " " /(^X/'^i) - /'(^s))) • 

The following theorem is proved analogously to Riemann Bilinear Identity 3.1. 
Theorem 6.1. For any two functions f,f: Q° — > C with periods 

A 1 ,...,A„B 1 ,...,B„A 1 ,...,A 9 ,B 1 ,...,B 9 and A[, . . . , A' g , B[, . . . , B' g , A[, . . . , A' g , B[, . . . , 1 
respectively, we have 

g 

(/, /') = £(A fc B' fe " BfcA * + AfcB 'fe " BfcA ^ 
fc=i 

One can prove analogously to First Existence and Uniqueness Theorem 2.3 that there is a 
unique (up to constant) Abelian integral 0g : Q° — > C of the 1st kind whose A-periods are given 
by the formula A^ = 5u, A& = Ski, where k — 1, . . . , g. Let B^, . . . , B^, B^, . . . , B^ be the B-periods 

of <f) l Q-. Q° — > C. The matrix (Ilg)^ := \{Q l k + B^) is the period matrix of the quad-surface Q. 

Problem 6.2. Generalize Convergence Theorem for Period Matrices 2.5 to quad-surfaces. 

Convergence of period matrices for quad-surfaces cannot be proved similarly to triangulated 
surfaces. The difficulty is the proof of a counterpart of Lemma 3.16 because for a quad-surface the 
natural "interpolation" I Q u is not continuous [25, Section 5.1]. 

6.2 Open problems 

Problem 6.3. Estimate the constants const-^., Const-^. in Convergence Theorem for Period Matri- 
ces 2.5 through the complex structure of 1Z only, e.g., through max{||rLft||, HIT^ 1 !!}. 

Problem 6.4. Generalize Convergence Theorem for Abelian Integrals 2.6 to Abelian integrals of 
the 2nd and 3rd kind. 

6.3 Underwater reefs 

Let us warn the reader on possible dangers arising on the way of proving convergence results. 

Remark 6.5. The proof of convergence of discrete period matrices suggested in [18, Proof of 3.1] 
is erroneous. It essentially relies on the assumption that the approximation of discrete analytic 
functions by discrete polynomials constructed in [ ] is uniform while the approximation is in fact 
pointwise and not uniform by [20, Section D]. 

Remark 6.6. A very different discretization of complex analysis is suggested in [33]. It is nonlocal, 
i.e., the notion of discrete analyticity depends on the values of a function at all the vertices, not 
just the neighbors. The proof of convergence of that discretization of period matrices suggested in 
[33, Proof of 7.2] is incomplete. It is stated there without any explanation that the convergence 
of discrete Abelian differentials implies the convergence of their periods. The norm in which 
the former holds is not specified. A careful analysis of the proof shows that the convergence of 
discrete Abelian differentials is proved in the L? norm, because the proof relies on the result [9, 
Theorem 2.8]. But the convergence in the L? norm is not enough to imply convergence of the 
periods (being integrals along certain 1-dimensional curves). 
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